Persistence Exponent for Diffusion: The Exact Solution 
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The persistence exponent 9o for diffusion equation cj>t = A(/!> with random gaussian initial condi- 
tions has been calculated for any dimension d. The value of 6q in the asymptotic limit for large time 
comes out to be d/4 . The result is at variance with the generally accepted values of 0.12,0.18 and 
0.23 for d = 1, 2, 3 respectively. 
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I. INTRODUCTION 

Persistence in non-equilibrium physics has been a field 
of active research for quite sometime now among both 
theoreticians and experimentahsts. Non-equilibrium 
processes may be of two types: stationary and non- 
stationary. Let us consider a stationary random pro- 
cess x{t) such that at t = this process is above a 
given threshold X = O(say). The interest lies in the 
probability P{t) for the process to have never crossed 
X upto time t — T. It is expected that for stationary 
processes the probabihty limT->.oo P{T) ~ exp[6'oT] . 9o 
is the persistence exponent where the subscript o in 
denotes the threshold 0. Infact, for non-stationary prob- 
lems, MmT^oo P{T) ^ T^^°. Non-stationary problems 
may be converted into stationary ones by a logarthmic 
substitution in time. 

Non-equilibrium aspects of statistical physics like diffu- 
sion[l,2], random walk[3], surface growth[4,5], have been 
studied in the context of persistence. Diffusion models in 
general have attracted a lot of attention in recent years 
in connection with their various features Apart from the 
simple diffusion equation which has been studied in the 
context of persistence, walks in 1-d random environments 
with or without bias [6] have been studied in details. 
Further, diffusion of a particle in one dimensional ran- 
dom potential with a small concentration of absorbers 
has been studied in the context of persistence [7]. There 
have been experiments too[8,9] to measure persistence 
exponents. 

One of the ways to calculate persistence exponent 
is through the Fokker-Planck equation-if it can be 
written down- subject to the required constraint on 
the probability. There are few exact calculations for 
the persistence exponent in the literature. The case 
of a simple random walk in one dimension gives the 
exponent 9o = -k- Even calculation of persistence expo- 
nents for gaussian processes may not be straight forward. 



We revisit the problem of simple diffusion, strongly 
non-Markovian in nature, given by the equation 0( = Acjj 
with random gaussian initial conditions. It has been 
an unsolved problem in statistical physics when for 
the first time results were reported by Majumdar 
et al[l] and Derrida et al[2]. Still, the problem of 
diffusion may require a better understanding in the 
context of persistence. The article tries to find an exact 
solution to the problem taking into account the fact 
that there is randomness only in the initial condition 
which propagates in accordance with the equation of 
motion. We conclude with a critical analysis of [1] and [2] 
in the context of the result obtained in the present article 



II. SIMPLE DIFFUSION 

The solution of the diffusion equation (j)t — A0 is 
uniquely determined by the initial condition. In the 
present problem, the initial condition is not fixed but 
is choosen from a distribution. The initial value of (p at 
every coordinate is choosen from a gaussian distribution 
and the initial values of (j) at any two coordinates are 
statistically independent 

In order to calculate persistence exponent we start with 
a positive at f = at the spatial coordinate Xo irre- 
spective of the values of (f> at other spatial coordinates. 
In the present case, the interest is in the probability Po(t) 
that (j) at a:o(say) remains positive without hipping even 
once after asymptotically large time t. The probability is 
expected to follow a relation Po{t) ^ t~^°, 6o being the 
persistence exponent. The random initial conditions for 
(j) are gaussian with 



{^{x,0))^0 



(1-a) 



(</)(xi , 0)ct,(x2 , 0)) = kS^"^ (Xi - X2) (1-b) 

where k is the variance. The solution for the diffusion 
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(f){x,t) 



d'^x'G{x-x',t)(P{x',0) 



(2) 



where G{x,t) = (47rt)~ ' exp(— x^/4t). We employ se- 
lective averaging over the initial field 0(x,O), where the 
averaging is done over i/i at all x at t = except at x — Xq. 
In other words, the averaging is done over all the initial 
configurations where at x = Xois kept fixed at a(say). 
Under the circumstances, the selective distribution, de- 
noted by s, is characterised by the moments. 
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(3-a) 



{(f>{xi,0)cf>{x2,0)),^ = {k+[a^-k]d^''\xi-Xo)}S^''\xi-X2) 

(3-b) 
It may be verified from (3 — &) that if xi = X2 7^ Xo, we get 
(1 — 6) and for xi = X2 = Xo, (3 — 6) gives a^ as expected. 
Using (2) and (3), we can calculate the moments of the 
random varible (f>{xo,t), 
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While evaluating the second order moment, we have used 
the relation in (3 — 6) . Hence the mean and the variance 
of the distribution for (j){xo,t), represented by /i and a^ 
respectively, are 



fx^i^ixo^ty^^iAnty/'c 
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(6-a) 



= k{ATTy'^2^'^/^-^^KdT{d/2)t-'^/'^ - kiAntyie-h) 

In the above equation kd denotes the angular integra- 
tion in d dimensional space while F represents the usual 
Gamma function. It may be mentioned that (t>{x, t) is 
gaussian irrespective of whether 4>{x' , 0) , the initial gaus- 
sian field, is correlated or not. This may be proved using 
Wick's theorem. In the present case, though, the initial 
field is uncorrelated. The S function distribution being 
the limitinp' case of a. o'a.iissia.n distribution, fh(T.„^ t\ is also 



a random gaussian variable. The expression for the con- 
ditional probability for starting at a and being between 
/3 and /3 + d/3 at time ti is 



P(/3|a)d/3 
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where /i = fi{a,ti) and cr — (y(t\). This probability con- 
siders all the paths that start from a to be between /? and 
/3 + dl3 after time ti including ones that flip enroute l3 as 
depicted in fig 1. Fig 1 is the projection of the trajectory 
of the system in the infinite dimensional $ — i space on 
to the 4>{xo) — t plane. 



4>(Xo) 




FIG. 1: 

A{0,a) represents the starting point and B{ti,f3), the 
destination. AB represents a path along which (j){xo) 
does not fiip and ADB is a typical path along which 
(j){xo) flips. Such paths have to be excluded. The prob- 
ability of reaching from A to the neighbourhood B at 
asymptotically large time ti without flipping is given by. 



P^ 



\a)dl3 ^ P{l3\a)dl3~ P{l3\~a)dl3 (8) 



The second term represents the probability of paths such 
as ADB originating from A(0, —a) and terminating in 
the neighbourhood of B at ti. Eqn (8), though has been 
used previously[10] , has a very different justification in 
the present case and holds good asymptotically. In the 
subsequent analysis we will consider a d dimensional lat- 
tice - lattice spacing being infinitesimally small- instead 
of continuum for the sake of notational convenience only. 
The reason for (8) follows. 

i) If an initial configuration at A of figl of the manuscript, 
given by Xi = {...ai, a, a2, ...}, is considered which 
takes (l){xo) to B , it may be inferred from (2) that 

f! + (4TTt)yya-^ -^ 



Xi ^ {.../ai,-a, /a2,...} (/ 
figuration at A which takes q 



/3-(47rt) "^'^a' 

{xo) to B. Hence there 
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and those from A to B . But it will be assumed that 
there is a one to one correspondence between the paths 
from A that flip to those from A to B.This may be used 
in the case of asymptotic t as a controlled approxima- 
tion since it improves with increasing t. In order to see 
this point let us consider a point B" (ii",/3)(not shown 
in the Figl) where ti" > ii. If Xf = {...ai,a,a2, ■■■} 
be the initial configuration corresponding to ADB , a 
path that flips, then X/' ~ {...fiai,a, fia2, ■■■} (/i = 
^tj_-|d/2 P-( TTti )^^^ a ^ jg ^j^g corresponding configura- 
tion from A to B" . The exact expression for /i con- 
tains a coordinate dependent term whose leading order 
behaviour for large t is 1. Since ii" > ti, we have /i > 1 
for sufficiently large ti. As a result, if the time coordinate 
at D be t^, (j>(xo,tD) < for the configuration X/' as a 
result of (2). Hence, one can conclude that the path cor- 
responding to Xf must have flipped at an earlier time 
than Id- Therefore, if a path from A ^^ B flips, the cor- 
responding path from A — > B" flips at an earlier time. 
Since ti" > ti, the 'number' of paths flipping while going 
from A -^ B" is more than those from A — > _B. Thus the 
'number' of paths from A ^ B that flip is a fraction F 
of those from A -^ B where F = 1 — 0{ti^'^) for large 
ti,a being some positive number. 

ii) Let 5i = {...,a, ...} be an initial configuration cor- 
responding to a path that flips before reaching B, a be- 
ing the value of (j){xo) at t = , the dots representing 
the values of (/> at all the other coordinates at t = 0. 
The value of only (j){xo) is now flipped so that we have 
S2 = {..., —a, •■•}■ The path corresponding to 6*2, AB' in 
fig 1 doesn't reach B but say B' at ti. In fig 1, the differ- 
ence between B and B' along the (^{xo) axis is 2a apart 
from a numerical factor depending on ti. As a matter 
of fact, in the asymptotic limit BB' — > 0. Initial (j) at 
every coordinate except at Xg may be suitably changed 
to bring B' back to B. Let the initial configuration thus 
obtained from 6*2 be denoted by 5*2. It may be easily seen 
that the paths corresponding to Si and 6*2 approach the 
same probability density in the asymptotic limit as / — )■ 1 
(see(i)) since the probability associated with any path is 
that of the corresponding initial configuration. As a re- 
sult, though 5*1 and ^2 do not have the exact probability 
density, we may multiply the probability of ^2 by a cor- 
rection term Tcor whose leading order term being 1 apart 
from higher order terms in ii . The value of 4> at Xo in S2 
or ^2 has been choosen as —a for the probability density 
at a and -a is the same on account of the distribution 
(1) being centered around zero. 

From (i) and [ii) we may say that that the probability 
of paths that flip while reaching B is P(/3| — ajFTcor- In 
the asymptotic limit (8) follows. The probability of not 
flipping after large time ti 






(9) 



value of t. Under the circumstances the second term on 

2 

the R.H.S of (6 — b) can be neglected. Further ^ ^ 
4>iO)t-'^/^ . Hence for 0(0) < f^/^, ^i^/cr^ < 1. The ex- 
pression (9) is evaluated using the identity[ll] 



dxexr>( 

^ 4/3 



7x) = y^cxp{/3j^)[l - erf(j^)] 



Hence we obtain 
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(11) 



In arriving at the above result the asymptotic expansion 
of 'error function' erf has been used for small argument. 
The expression for Po{t) is 



pQO 

Po{t)= daP+{a)^^.Q{a) 



(12) 



We -would no-w calculate (Q^ for asvmntoticallv larce 



where Q{a) is the gaussian distribution for a with 
variance k as mentioned at the beginning. Ii k « f^^^, 
it may be inferred from (11) and (12) that Po{t) ^ t^'^/^. 
This gives 60 = d/4 

The answer of d/A is at variance with that of [1] and 
[2]. The difference may be accounted for in their appli- 
cation of the 'independent interval approximation' (HA), 
a Markovian one. 



III. DISCUSSION 

The persistence exponents for the problem of simple 
diffusion with random, gaussian initial conditions were 
evaluated in [1] and [2] to be Oo = 0.12,0.18,0.23 for 
d ~ 1,2,3. The central feature in the approach of [1] 
and [2] is the approximation known as the 'independent 
interval approximation' (HA) and the subsequent results 
are an artifact of HA. The results, as obtained using 
the approximation, have been confirmed, as claimed 
by [1,2], by simulation results. It is argued that both 
the approximation(IIA), an uncontrolled Markovian 
approximation, as well as the algorithm (Euler) for the 
simulation are inappropriate for the present problem 
and hence the validity of the results may be questioned. 
The problem is non-Markovian and the Markovian 
approximation (HA) ignores the fundamental feature 
of the problem. Regarding the simulation, the use 
of Euler's algorithm , according to literature [8], is 
untenable for diffusion on a large lattice because of 
certain serious complications which are very specific 
to the diffusion equation. The non-Markovian nature 
as well the complications surrounding the simulation 
of the present problem have been elaborated upon below. 

i) A Markov process has no memory. In the present 
nrnhlem it is encoded in the initial conditions (h(x.^ fX\ 



through <j>(x,t) = J d'^x'G{x — x' ,t)(f>{x' ,0), where 
G{x — x',t) is the kernel for the diffusion equation, 
as to at what time at a particular coordinate will 
flip making the process strongly non-Markovian. The 
problem of simple diffusion is not only non-Markovian 
but is different from other problems encountered in 
the field of persistence for only the initial condition is 
stochastic and the subsequent evolution is deterministic. 
Hence, the problem needs to be addressed in a different 
way. The probability of not flipping after time t is the 
sum over the probabilities of the initial configurations 
corresponding to paths that do not cross zero (at 
Xo) even after time t. The approach in the present 
article depends crucially on the initial condition and it's 
distribution through (7) and (8). 

ii) The diffusion equation is unique in the sense that the 
Euler's method cannot be used whenever the lattice size 
becomes large (10^). The stability criterion in the case 
of Euler's equation for unit diffusion constant is given 



by 



2At 



< f [12] where Ai is the time step and Ax is 



the lattice constant. If Ax — 1, then At < i. Hence the 
time step is invariably small in this method. It is also 
well known that T ^ L^ , T is the time scale of diffusion 
and L the system size. If the system size is 10®, the 
corresponding time scale comes out to be 10^^. Hence 
in order to probe the region of interest where things 
happen, we must have t > 10^^. Hence t = 10^ cannot 
be considered asymptotic, as has been considered in 
[1,2], This is where Euler's method fails for At is small 
as mentioned before. It may be concluded that if Euler's 
method is used, the computation of persistence exponent 
becomes inconceivable for it is not possible to cross a 



temporal region of width ^ 10""^^ using small Ai. Crank- 
Nicholson(CN) or Backward Timc(BT) methods [12] are 
superior to Euler for there is no constraint on At . One 
can use use the Fourier space version of the diffusion 
equation for the simulation. Further, as claimed in [1], 
the result for one dimension has been arrived at after 
averaging over 17 runs. Infact, the number is very low 
for obtaining convergence. Rather, averaging over such 
small number of runs is marked by large fiuctuations 
about the mean. 



Simulation has been performed by us for the problem 
using spectral method with d — I, diffusion constant 
D — .001 for a lattice of size 10"* and 800 random ini- 
tial values to observe that till time 10^ one may obtain 
6*0 = 0.12 but between 10"^ and 10^ there is no flipping 
i.e 6*0 = though there are paths, 120 out of 800, that 
have not flipped. Hence, till the point we have done the 
simulation, it tells us that the exponent is and not 0.12 
as is generally claimed. This may hint towards a cross- 
over which would be clear if the simulation is run for a 
long enough time(t ^ 10^^) and a much larger sample size 
which we didn't. It took us one month to run the parallel 
code using 40 processors on a cluster . It may be added 
here that there are other numerical results [13] which ap- 
pear to validate [1,2]. In [13] the averaging has been done 
such that the fact of stochasticity being present only at 
t ~ 0, the central aspect of the problem, gets ignored. 
Though there is some literature connected to the prob- 
lem of persistence in diffusion, it appears that a few issues 
still need to be addressed. 



[1] Satya N. Majumdar, Clement Sire, Alan J. Bray, and 

Stephen J. Cornell, Phys. Rev. Lett. 77, 2867(1996). 
[2] B. Derrida, V. Hakim, and R. Zeitak,Phys. Rev. Lett. 

77, 2871(1996) 
[3] J.M. Schwarz and Ron Maimon, Phys. Rev. E. 64,016120 

(2001) 
[4] J.Krug, H.Kallabis, S.N.Majumdar, S.Cornell, A.J. Bray, 

and C.Sire, Phys.Rev.E.56, 2792(1997) 
[5] M. Const ant in, S.Das Sarma, C.Dasgupta, 

Phys.Rev.E.69, 051603(2004) 
[6] Daniel S. Fisher, Pierre Le Doussal and Cecile Monthus, 

Phys. Rev. E 59,4975(1999) 
[7] Pierre Le Doussal, arXiv: condmat/0906.0267 



[8] B.Yurke et al., Phys.Rev.E56,R40(1997) 

[9] W.Y.Tam et al., Phys. Rev. Lett. 78, 1588(1997) 

[10] S.Chandrasekhar, Selected Papers Vol 3: Stochastic, Sta- 
tistical and Hydromagnetic Problems in Physics and As- 
tronomy, The University of Chicago Press, 1989 

[11] I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Se- 
ries and Products, Academic Press, 1980 

[12] W.H. Press,S.A. TeukolskyW.T. Vallering and B.P. 
Flannery, Numerical Recipes in C, Cambridge University 
Press, 1992 

[13] George C.M.A. Ehrhardt and Alan J. Bray, Phys. Rev. 
Lett. 88, 070601(2002) 



